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In addition to density perturbations, inflationary models of the early universe generally predict 
a stochastic background of gravitational waves or tensor fluctuations. By making use of the infla- 
tionary flow approach for single field models and fitting the models with Monte-Carlo techniques 
to cosmic microwave background (CMB) data from the Wilkinson Microwave Anisotropy Probe 
(WMAP), we discuss the expected properties of the gravitational wave background from inflation at 
scales corresponding to direct detection experiments with laser interferometers in space. We comple- 
ment the Monte-Carlo numerical calculations by including predictions expected under several classes 
of analytical inflationary models. We find that an improved version of Big Bang Observer (BBO- 
grand) can be used to detect a gravitational wave background at 0.1 Hz with a corresponding CMB 
tensor-to-scalar ratio above I0 -4 . Even if the CMB tensor-to-scalar ratio were to be above 10 -2 , we 
suggest that BBO-grand will be useful to study inflationary models as the standard version of BBO, 
with a sensitivity to a stochastic gravitational wave background Qcwh 2 > 10~ lr , will only allow a 
marginal detection of the amplitude while leaving the tensor spectral index at 0.1 Hz unconstrained. 
Also, inflationary models with a large tensor-to-scalar ratio predict a substantial negative tensor 
spectral index such that the gravitational wave amplitude is damped at direct detection frequencies. 
We also discuss the extent to which CMB measurements can be used to predict the gravitational 
wave background amplitude in a direct detection experiment and how any measurement of the am- 
plitude and the spectral tilt of the gravitational wave background at direct detection frequencies 
together with the CMB tensor-to-scalar ratio can be used to establish slow-roll inflation. 

PACS numbers: 98.80.Bp,98.80.Cq,04.30.Db,04.80.Nn 



I. INTRODUCTION 

The data from high precision cosmological experiments 
relating cosmic microwave background (CMB) radia- 
tion [l| and large-scale structure mass distribution from 
galaxy redshift surveys have clearly demonstrated that 
the inflationary paradigm is the most preferred de- 
scription of the origin of density perturbations. While 
expectations are generally consistent with existing ob- 
servations, we still lack a complete understanding of un- 
derlying physics related to inflation including the shape 
and the amplitude of the potential responsible for super- 
luminal expansion. 

The cosmological observations, so far, have been re- 
stricted to studying the density, or scalar, fluctuations 
in the universe both through CMB and large-scale struc- 
ture. In addition to scalar perturbations, tensor fluc- 
tuations in the form of a primordial gravitational wave 
background are expected to be generated during infla- 
tion [3j. If detected, the amplitude of the gravitational 
wave background provides a direct estimate of the poten- 
tial height of the inflaton field when relevant modes exit 
the horizon. This background is now the focus of a large 
number of ongoing and upcoming CMB polarization ex- 
periments since gravitational waves leave a distinct sig- 
nature in the polarization pattern of CMB photons in 



the form of B- or curl-modes Q. The Wilkinson Mi- 
crowave Anisotropy Probe (WMAP; data have now 
constrained the the tensor-to-scalar ratio, tcmb, denot- 
ing the tensor power spectrum amplitude divided by the 
amplitude of the scalar perturbation power spectrum, to 
be below 0.38 at CMB scales (at 95% confidence level 
i). 

The planned next generation NASA mission for high 
precision polarization observations ( CMBpol or Inflation 
Probe) has now set a detection goal for tcmb above 
10 -2 In the future, with an ambitious experiment, 
this limit could be further improved, though confusion 
from known polarized foregrounds involving dust and 
synchrotron will limit CMB studies to a level above a 
few times 10 -4 While all-sky maps of polarized syn- 
chrotron emission are now available with low-frequency 
channels of WMAP, the aforementioned uncertainty is 
due to the lack of all-sky maps at frequencies dominated 
by polarized dust. With data from Planck, it is very 
likely that we will have a more exact estimate on the 
limit to which a tensor background can be probed with 
next generation CMB polarization experiments. It is now 
clear that the limiting factor for CMB observations is 
not the ultimate lensing confusion that was previously 
discussed in the literature with a minimum detectable 
tensor-to-scalar ratio around 10 -5 KM. 
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Beyond CMB, there is now a growing interest to di- 
rectly detect the relic gravitational wave background 
from inflation with laser interferometers in space with 
test masses separated by a few thousand to 50,000 kilo- 
meters. Plans for such studies include NASA's Big Bang 
Observer (BBO) ^} and the DECI-hertz Interferome- 
ter Gravitational Wave Observatory (DECIGO) in Japan 
[Tl| . Other studies that explore possibilities for future 
gravitational wave observations beyond first generation 
Laser Interferometer Space Antenna can be found in 
Ref. 01 • Based on the expected foreground confusion 
and technological improvements, current concept studies 
aim the frequency regime between 0.1 Hz to a few Hz 
[lfij . Here we will use 0.1 Hz with a corresponding scale 
k ~ 6.47 x 10 13 Mpc -1 as the fiducial frequency for direct 
detection experiments. As CMB measurements are sen- 
sitive to density and tensor perturbations at horizon-size 
scales, at k ~ 0.002 Mpc -1 , the two detection techniques 
are roughly sensitive to a gravitational wave background 
at wavelengths separated by 16 orders of magnitude. 

As discussed in a variety of recent papers 0, 0, 0, 
lili IT?! ] , the large lever arm produced by combining infor- 
mation at CMB and direct detection scales is expected 
to pin down inflationary models better than data from 
any one of the two scales alone. Recent studies have 
discussed the expected signal for direct detections either 
through exact calculations using specific analytic models 
[l5j | or through the use of numerically generated Monte- 
Carlo models 0, under the inflationary flow equa- 
tions pbl l2l| . While general predictions exist and the 
expected level of the background is understood in terms 
of the CMB tensor-to-scalar ratio, past studies have paid 
little attention to issues of connecting CMB information 
with the gravitational wave background data at direct de- 
tection frequencies to improve inflationary models. While 
the long lever arm connecting the two scales has been 
advocated to help discern inflationary model space, as 
we discuss here, such a large lever arm also complicates 
a simple analysis based on a power-law extrapolation 
of tensor and scalar power spectra from CMB scales at 
k = 0.002 Mpc -1 to direct detection at k = 6.47 x 10 13 
Mpc -1 . 

Furthermore, while past studies have discussed the 
possibility for a direct detection of the inflationary back- 
ground with future laser interferometers [2(|, the sen- 
sitivity requirements for a direct detection experiment 
have not been fully characterized. We will try to de- 
termine how sensitive direct detection experiments will 
need to be and put these sensitivities in context with 
CMB polarization measurements. We also extend previ- 
ous analysis that only concentrated on the gravitational 
wave background amplitude [l^ to also discuss the de- 
tectability of the tensor spectral index at these frequen- 
cies. Our discussion is similar in spirit to the recent study 
in Ref. |16| . but we improve the analysis by considering 
both numerical as well as analytical 15] models of infla- 
tion, combining separate analyses in the literature into a 
single discussion. 



The combination of the tensor amplitude and the tilt 
can be used to validate inflationary models as these ob- 
servables are expected to satisfy certain consistency re- 
lations. Similar tests are discussed at CMB scales in the 
literature, though these relations are poorly constrained 
due to the fact that tensor tilt is generally hard to mea- 
sure at CMB scales j^. As an application of our models, 
we extend the discussion in Smith et al. |l6j | to consider 
the possibility to test the inflationary consistency rela- 
tion. This is done in Smith et al. under the assumption 
that one can use the tensor tilt from direct detection ex- 
periments as a proxy for the tensor tilt at CMB scales. 
Using numerical models, we show that for general infla- 
ton potentials the tensor tilt at k = 6.47 x 10 13 Mpc -1 
differs from the tensor tilt at k = 0.002 Mpc -1 by more 
than a factor of ten. Furthermore, beyond simply using 
a direct detection experiment to confirm a detection of 
tensor modes at CMB scales and thereby improve infla- 
tionary model selection, we also explore the usefulness 
of a highly sensitive direct detection experiment that is 
also capable of measuring the tensor tilt at direct detec- 
tion frequencies. Unfortunately, while experiments such 
as ultimate-DECIGO can reach gravitational wave back- 
ground amplitudes with the corresponding CMB tensor- 
to-scalar ratio of 10 -6 and above, tensor tilt will remain 
unconstrained for most of the models with CMB tensor- 
to-scalar ratios below 10~ 3 and tilts with magnitude such 
that \n T \ < 0.01. 

This paper is organized as follows: In the next Sec- 
tion, we outline our calculations that we use to predict 
properties of the gravitational wave background from in- 
flationary models. We consider two separate approaches, 
one involving Monte-Carlo models of the inflationary flow 
equations, and the second using analytical families of in- 
flationary models that are consistent with observations 
at CMB scales. In Section III, we discuss our results and 
comment on the sensitivity requirements for a direct de- 
tection experiment to reach similar sensitivities as those 
expected with CMB polarization observations. We con- 
clude with a summary of our main results in Section IV. 



II. CALCULATIONAL METHOD 

A. Monte-Carlo Reconstruction of Inflationary 
Flows 

To model the evolution of single field inflationary mod- 
els, we study dynamics through the hierarchy of flow 
equations involvin g th e generalized "Hubble Slow Roll" 
(HSR) parameters [2(J, l2lL |22| . The advantage of the in- 
flationary flow method is that one can study the generic 
behavior of the inflaton field without assuming a par- 
ticular shape for the potential. We outline the main 
equations of this approach here and refer the reader to 
Refs. jHlll for further details. 

To begin, we note that the motion of the scalar field in 
a cosmological background is given by (f> + 3H<f> + V' (4>) = 
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where an overdot corresponds to the time derivative, 
a prime denotes the derivative with respect to <j) and 
H = (a J a) is the Hubble parameter. Assuming the infla- 
ton field dominates the energy density, the Einstein field 
equations are 
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As discussed in Refs. 
be written as 
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p3l |. the equations of motion can 



4tt 
12tt 



(2) 



M pl z M p i 

where the Hubble parameter is now a function of the 
field 4> instead of time under the assumption that (j> varies 
monotonically with time. 

As usual, dynamics are described in terms of the HSR 
parameters that are defined in terms of the derivatives of 
H with respect to cj> [24|: 
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and for c omp arison with similar calculations in the liter- 
ature [13, |2(J, Wi H3 we identify the slow roll parameters 
1 \h and 2 Xh as r\ and £, respectively. We also define 
cr = 2?7 - 4e. 

The trajectories of these HSR parameters are governed 
by a set of coupled first order differential equations [2£j, 
l2ll ] written in terms of the number of e-folds, N, before 
the end of inflation, with the convention that N increases 
as one goes further back in time, as 
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X H ) + (^Xh) (£ > 1). 



As written, the evolution of a given slow roll parame- 
ter depends on a slow roll parameter that is one order 
higher. While the coupled differential equations involve 
an infinite hierarchy, in practice, this hierarchy must be 
truncated at some order M such that M+1 A# = 0. From 
Eq.EJ with d^ M+ ^H/d(j) {M+ ^ = 0, one can then write 
H ((f) as a polynomial of order M + 1 with |2^| 
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Further, from the definition of e(0), one can also write 
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when the coefficients Ai, with i > 1, are related to initial 
values of HSR parameters 
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with A\ = y47reo and the sign of A\ determining the 
direction the field is rolling. 

Once e(N) is determined as a function of the number of 
e-folds before the end of inflation N, the relation between 
number of e-folds and the field 6 is determined from 



d(j> 
dN 



pi 



while the Hubble parameter is established from 

1 dH _ 
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With e(4>) and H((f) reconstructed, one can also obtain 
the inflaton potential through the Hamilton- Jacobi equa- 
tion 



(10) 



Using equation (10) and the expressions in (Q, one can 
write (a/a) = H 2 ((j))[l — e((j>)). Inflation is defined as 
(d/a) > 0, which happens as long as e((j)) < 1. We use 
this criteria in each of the numerically generated models, 
such that the end of inflation occurs when e(</>) — > 1. 

Note that there is an overall constant of Hq here that 
is not set by the differential equations. This constant can 
be established based on the observed amplitude of cither 
tensor or scalar fluctuations at a specific scale. Here, 
we will normalize our models based on the amplitude 
of density fluctuations at scales corresponding to CMB 
observations as determined by WMAP data. Since mod- 
els are normalized to CMB, we pick the fiducial physical 
scale that corresponds to <Pcmb — to match the usual 
pivot point of CMB observables at fccMB = 0.002 Mpc -1 . 
Beyond this pivot point, the physical wavenumber is as- 
sociated with a value of </> through 
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Our convention is such that cf> > corresponds to scales 
smaller than kcMB (k > 0.002 Mpc -1 ) such that one is 
going ahead in time. 
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At a given wavenumber, the observable power spectra 
of tensor and scalar perturbations are given as [34j 
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when x = -2 + In 2 + 7 and 7 = 0.5772156649 is the 
Euler-Mascheroni Constant. This allows us to set the 
overall normalization Hq such that at fccMB we take 



P s (k = 0.002) = (2.45 ± 0.23) x 10" 



(13) 



consistent with WMAP p] . This normalization is equiv- 
alent to setting: 
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and for most practical purposes the term within the 
square bracket can be ignored [3|| . With our choice that 
CMB scale is at = 0, in equation (5), H = Hqmb- 

Instead of amplitude at power-spectra at each k, the 
observables are generally described in terms of the power- 
law variables given as 
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where C = 4(ln2 +7) - 5 = 0.08145. For simplicity, we 
also discuss the tensor-to-scalar ratio defined at a partic- 
ular scale as 



be used to validate inflationary models. We will discuss 
the accuracy to which this relation is valid at both CMB 
and direct detection scales. Though we can only expect 
a measurement of the tensor amplitude at 0.1 Hz with 
direct detection experiments, we will discuss the relation 
between tensor-to-scalar ratio at CMB scales tcmb and 
the same ratio at direct detection frequencies, ro.iHz- Due 
to the scale dependence of the scalar and tensor tilts, 
these two ratios are not expected to be equal for a general 
inflationary potential. 

To generate a family of models, we follow the usual 
recipes outlined in the literature 0, l2ll . We take 
the initial conditions to be sampled with in the parameter 
range of 

eo € [0,0.8] 

cr e [-0.5,0.5] 

fo e [-0.05,0.05] 

e X H \o G [-0.025 x 5~ e+3 , 0.025 x 5~ e+ % (3 <£< 10) 

11 ^H|0 = 0, 

such that the hierarchy is truncated at the tenth order 
in HSR parameters. We solve the flow equations for each 
set of numbers randomly drawn within the above ranges. 
Assuming each of the ranges is uniform, we establish 
[e(N), e Xh(N)]. We require that inflation last between 
46 and 67 e-folds, roughly consistent with known con- 
straints [30II3H . With e(N), we also solve for the relation 
between <f> and N and between <f> and k. 

The latter establishes primordial power spectra as a 
function of k. The relations also establish the power-law 
observables at a given wave number, though to compare 
with observations we restrict them to the values at CMB 
scale, with A: C mb = 0.002 Mpc" 1 , and fc dir = 6.47 x 10 13 
Mpc -1 . The large number of e-folds between these two 
scales, A7V = ln(fcdh //ccMB) ~ 38, is expected to improve 
constraints on the inflationary model since the potential 
will be probed at two largely different field values |15|. 

Using the above initial conditions, we ran close to 20 
million individual models |3jJ. Out of these 20 million 
models, roughly 20,000 models are found to fall within 
WMAP 95% confidence level constraints of tcmb < 0.38 
and 0.92 < n s < 1.06. The latter are suggested from 
recent cosmological data analyses (see e.g. 0) and 
broadly describe the allowed parameter ranges within the 
2er confidence level. In Fig. 1, we summarize the (n s ,r) 
plane and the distribution of model points. We will dis- 
cuss these and other results in Section III. 
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which in terms of slow-roll parameters is 

r= 16e [1 - C(cr + 2e)] . (17) 

To the first order in slow roll parameters, r w — 8n t , 
which is considered to be a consistency relation that can 



B. Analytical Models of Slow- Roll Inflation 

In addition to the numerical method to describe infla- 
tionary predictions, we also make use of several analytical 
models of slow-roll inflation to connect CMB observables 
with those at direct detection experiment levels. Unlike 
the previous subsection where flow equations allow us 
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FIG. 1: The region in (n a ,r) plane consistent with WMAP CMB data at the 2a confidence level with tcmb < 0.38 and 
0.92 < n s < 1.06, where rcMB is the tensor-to-scalar ratio and n s is the scalar spectral index, or tilt. Both these parameters are 
determined at CMB scales with k = 0.002 Mpc -1 . Inflationary models corresponding to each model point shown here also give 
rise to N > 46 e-folds in the scale factor before the end of inflation. For comparison, we also plot the parameter space occupied 
by four models of inflation outlined in Section IIB: power-law (solid blue line), chaotic (dotted magenta), symmetry breaking 
(dot-dashed cyan), and hybrid (dashed yellow) models. The two horizontal lines outline the expectation on the tensor-to-scalar 
ratio from CMB studies. The two values are chosen to the target goal of the next generation CMB polarization experiment 
(rcMB > 10~ 2 ) and an optimistic limit related to polarized foreground confusion (ccmb > 10~ 4 ). A more realistic estimate on 
the foreground confusion will suggest that CMB studies will be limited anywhere between a tensor-to-scalar ratio of 10~ 3 to 
KT 4 . 



to describe the dynamics without specifying a particular 
shape for V((j>), with analytical models, we are required 
to specify the potential explicitly. For easy comparison, 
we take the same four families of models as those consid- 
ered in Ref. . These models describe both single field 
models and hybrid inflation models. The four potentials 



considered here are 

V{4>) = V e~ p ' t ' /Mpl , Power - law 

Chaotic 
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, Hybrid (18) 
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FIG. 2: The running of the scalar spectral index against the spectral index at CMB scales. We did not constrain the running, 
although a majority of the models fall below 10 -2 and have positive running. The horizontal line at = 0.007 shows the 
WMAP 2a upper-limit on the running of the scalar spectral index. 



and describe power-law, chaotic, symmetric breaking, 
and hybrid inflation, respectively from top to bottom. 
The calculation involves simply matching the observables 
at CMB scales (n s , r) to free parameters in each of the 
potentials, such as a and v, and also uses the number of e- 
folds between CMB horizon exit and the end of inflation. 
The amplitude of density perturbations Ps (fccMB) fixes 
the overall potential normalization given by Vq. With 
Vq and free parameters fixed, and taking AN = 38, we 
calculate the observables in terms of f^cw^ 2 (see below) 
and n-T and, in cases where running is not zero, ax at 
the field value <fi corresponding to direct detection exper- 
iments. We do not reproduce the analytical expressions 
associated with the calculation of predictions under these 
models as these are either readily available in Ref. [15( or 



derived through equations outlined there. 



C. Direct Detection Observables 

While CMB measurements involve Ps{k) and Pr(k), 
direct detection measurements involve Pt (k) only. As we 
discuss later, the lack of information on the primordial 
scalar spectrum at direct detection frequencies compli- 
cates simple tests and studies related to slow-roll inflation 
using data from direct detection experiments alone. Fur- 
thermore, since the observations involve the relic back- 
ground today, one must also account for the evolution 
of the tensor mode from high redshifts to low redshifts. 
This evolution of a single gravitational wave mode can 
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FIG. 3: The set of potentials that give rise to points shown in Figure 1 from the Monte-Carlo inflationary flow analysis. Here 
4> = corresponds to observations at CMB scales with k = 0.002 Mpc -1 , while <f> > corresponds to scales smaller than CMB. 
The potentials are subdivided based tensor-to-scalar ratio at CMB. Potentials that give rise to large tensor-to-scalar ratios with 
fCMB > 10~ 2 cannot simply be described by simple power-laws or any particular shape. 



be written through the massless Klein-Gordon equation 

d 2 hi. 1 da dhu , o, . . 

-TT + 2-- r -7 1 + k 2 h k = 19 
dr z a dr dr 

with the boundary conditions /ife(0) — Pr(k) 1 / 2 and 
^fe(O) = 0. Here, r is conformal time. We define 
g'fc(r) = /ifc(r)a(r). With this definition, we are able 
to rewrite the above equation, 

d 2 g k ( , 1 d 2 a\ , x 

-d^ + { k -a^> = ' (20) 

There are two limiting behaviors for g k : before horizon 
entry g k oc a — > h k — constant; after horizon entry g k will 
oscillate which implies h k will oscillate with an amplitude 



decreasing as 1/a |19| . Therefore, the current spectrum 
of gravitational waves is determined by the primordial 
power spectrum and by the rate at which scales enter 
the horizon (i.e. the evolution of the scale factor): h — 
hodk/a. During radiation domination, H oc a~ 2 , so that 
k = l/dk and during matter domination H oc a~ 3 / 2 , so 

1 /2 

that k = l/a k . From these relations, we find 

hk oc k^ 1 (Radiation Domination) (21) 
hk oc k~ 2 (Matter Domination) . (22) 

These scalings are generally quantified in terms of the 
transfer function Tx(k) for tensor modes 19], such that 
the tensor power spectrum is related to the one generated 
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FIG. 4: The relation between the number of e-folds N and <p, the value of the inflaton field, for representative sampling of 
potentials shown in Figure 3. In the left panel, we make use of potentials that give rise to t"cmb > 10 -2 , while in the right panel, 
we concentrate on the relation between N and (j> for potentials that lead to tomb < 10~ 4 . In the case of potentials that lead 
to large ?"cmb , V{tf>)'B are shaped such that there is a significant, and an abrupt, increase in the number of e-folds over a small 
variation in <f>. The <f) values related to this abrupt changes in N correspond to the field values where the inflaton potentials in 
Figure 3 flatten. On the other hand, N versus d> relations for potentials that results in a small ?"cmb are gradual. A behavior 
similar to this was already described in Ref. [17| based on the relation between H(N) versus N such that models that lead to 
large tensor amplitude have abrupt changes in H(N). The horizontal lines indicate the number of e- foldings related to CMB 
(top two lines) and direct detection (bottom two lines) observations. 



during inflation as 

P T (k,t = t ) = P T (k,t = t oo )T*(k). (23) 

We define the ratio of energy density of the gravitational 
wave background to the critical closure density as 



^Gw(fc) 



1 dp GW 
p c dink 



(24) 



and write 



n GW (k)h 2 = 



m 2 



\h k \ 2 )=A GW P T (k), (25) 



where A GW = 2.7 '4 x 10~ 6 at f = 0.1 Hz |15| and comes 
from taking a frequency average of the tensor transfer 
function Ty(fc). The mapping amplitude ^4gw ignores 
the anisotropic stresses in the cosmic fluid, among which 
free-streaming of relic neutrinos is important 29] . The 
additional damping of tensor modes, however, is only re- 
stricted to modes that enter the horizon after neutrino 
free-streaming or the gravitational wave background at 
frequencies below 10 -11 Hz. Thus, our calculations are 
not affected by ignoring anisotropic stresses. We do note 
that other exotic processes may damp gravitational waves 



at frequencies of 0.1 Hz 28], but the physical processes 
are mostly speculative and we do not consider them to 
be important for this calculation. 



In addition to the amplitude of the gravitational wave 
background captured by f^cw^ 2 , recent analyses suggest 
that it may also be possible to measure the tilt ut of 
the gravitational wave background at direct detection fre- 
quencies [2(|. Thus, we discuss the observables in terms 
of both ficw^ 2 and Ut here. For direct detection ex- 
periments, we consider three possibilities here involving 
BBO and an improved version of BBO (BBO-grand) as 
well as DECIGO. These experiments and their sensitiv- 
ities are summarized in Ref. pp| . and to be consistent 
with those calculations, we make use of the same quoted 
sensitivities here. We do not consider a version of BBO 
called BBO-lite with a sensitivity to gravitational waves 
only down to fl G yjh 2 of 10~ 15 . As we find later, such an 
experiment is unlikely to result in any detection of the 
gravitational wave background as current limits already 
rule out a background just above this level. 
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FIG. 5: The amplitude of the stochastic gravitational wave background, Qcwh 2 , versus the gravitational wave frequency. The 
curves, subdivided by the tensor-to-scalar ratio at CMB scales, represent the expectations based on inflationary models selected 
from solving inflationary flow equations and consistent with WMAP+SDSS data. The models that lead to large tensor-to- 
scalar ratios generally have large spectral indices for the gravitational wave background resulting in a damping of the power 
at frequencies related to direct detection experiments. We also show the detector sensitivities for a 5 year-long integration of 
some of the options discussed in the literature. These experiments, from top to bottom, are BBO, an improved version of BBO 
(BBO-Grand), and an extremely optimistic version of the DECIGO (Ultimate DECIGO). 



III. RESULTS 



A. CMB Predictions and Comparison 



Figure 1 summarizes the plane of tensor-to-scalar ra- 
tio versus scalar tilt at CMB scales (k = 0.002 Mpc" 1 ), 
where points represent results of the Monte-Carlo pro- 
cess using the flow equations. We constrain the parame- 
ter space to be tcmb < 0.38 and 0.92 < n s < 1.06 at the 
95% confidence level following the results of Ref. [l| and 
Ref. @. 



In addition to lying within this range, we also require 
that each of the selected models generates between 46 and 
67 e-folds before the end of inflation. The two horizontal 
lines represent the expected limits from CMB observa- 
tions: the published goal of the next generation NASA 
CMB polarization mission with a tensor-to-scalar ratio of 
10 -2 0| and an optimistic limit on the tensor-to-scalar 
ratio based on polarized foregrounds (10~ 4 ). This limit 
is somewhat unknown and estimates suggest that it is 
between 10~ 3 and 10~ 4 ||, depending on the frequency 
selection, the sky area targeted for CMB observations, 
and the unknown polarized intensity of dust. Once high 
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FIG. 6: The amplitude of the gravitational wave background 
at 0.1 Hz flGwh 2 plotted in terms of the tensor-to-scalar ra- 
tio at CMB scales (k — 0.002 Mpc -1 ). For reference, we also 
show the minimum background amplitude that can be de- 
tected with a signal-to- noise ratio of unity with BBO, BBO- 
Grand, and Ultimate DECIGO. The relation between QGwh 2 
and rcMB is such that one can establish an upper limit on the 
expected amplitude of the gravitational wave background at 
a given tensor-to-scalar ratio. The dotted line shows the re- 
lation, which we discuss in equation 27. An experiment such 
as BBO-Grand can probe slow-roll inflationary models with 
tensor-to-scalar ratio greater than 10 -4 . 



frequency all-sky CMB polarization maps become avail- 
able with Planck, this foreground- limited tensor-to-scalar 
ratio will be better established. To guide our discussion, 
here, we take the lower limit of 10 -4 . 

In addition to inflationary flow equation predictions, 
we also show the expected range allowed by the model 
parameter space of several families of analytical models of 
inflation. These include the power-law, chaotic, symme- 
try breaking, and Hybrid inflation models. The allowed 
ranges we hig hlight here are consistent with previous cal- 
culations [15| and both analytical and Monte-Carlo nu- 
merical models are normalized to the same amplitude of 
the scalar perturbations at k = 0.002 Mpc -1 . As is well 
known, hybrid inflation models generally allow n s > 1, 
while other families generally fill up the parameter space 
with n s < 1. 

For comparison, in Figure 2, we plot the scalar spectral 
index against the running of the scalar spectral index, a s , 
at CMB scales. The top line indicating an arrow pointed 
to low values of a s is roughly the upper limit on the 
running of scalar tilt with a value of 0.007 at the 95% 
confidence level using WMAP+SDSS data 0. At the 



low end, the running is allowed to be as low as -0.13 and 
all our data points lie above this limit. While some of our 
models lie above a s — 0.007, we do not impose a selection 
based on a s , as we have seen from the literature that 
the determination of a s is largely subject to assumptions 
coming from model fitting. As an example, in Ref. |32| . 
a s is constrained to be below 0.01 at the la confidence 
level with a prior of 30 e-folds or more between CMB and 
the end of inflation. 

The corresponding potentials allowed by the inflation- 
ary flow analysis in the (n s ,r cmb) plane are shown in 
Figure 3. The potentials are subdivided based on the 
tensor-to-scalar ratio at CMB scales. In Figure 3 <f> = 
corresponds to CMB scales while <f> > represents scales 
below CMB. The potentials are truncated at the end of 
inflation when e reaches a value of unity, and again, all 
these potentials lead to between 46 and 67 e-folds of infla- 
tion. As shown, potentials that give rise to large tensor- 
to-scalar ratios with tqmb > 10~ 2 cannot be described 
by simple power-laws or with any particular analytical 
shape. In fact, almost all of the potentials cannot be 
described by simple analytical forms. 

In general, inflationary potentials with tcmb > 10~ 2 
have A(j) > M p \, where A<j) corresponds to the number of 
e-folds from CMB horizon exit to end of inflation. This 
result is well known in terms of the Lyth-bound |33J . In 
Fig. 4, we summarize the relation between the number 
of e-folds N and the field value cf> for potentials that give 
rise to tcmb > 10~ 2 (left panel) and tcmb < 10~ 4 (right 
panel). As can be seen in Fig. 3, potentials that lead to 
large tensor-to-scalar ratio at CMB scales have structure 
such that towards the end of inflation potentials begin 
to flatten followed by a drop in the amplitude as £ ^ 1. 
Since the number of e-folds between two field values can 
be written as 

AW)~-^r^>, (26) 
M p i J fa V'{<p) 

and since V'(cf>) becomes smaller towards the end of in- 
flation, one finds a large change in the number of e-folds. 
In Fig. 4 we show some example relations between N and 
(f> which highlight the fact that most potentials that lead 
to rcMB > 10~ 2 behave such that the last 20 to 30 e- 
folds of inflation happen suddenly over a small duration 
of A(j). The models with steep slopes around = 2 are 
prime examples of this behavior. Similar potential shapes 
have been noted before and explained in terms of the re- 
lation between H(N) and N in Ref. 01 ■ As discussed 
there, direct detection experiments that are only sensi- 
tive to gravitational wave backgrounds with large tcmb 
values will be probing potentials with this abrupt behav- 
ior. None of the analytical models have behaviors like 
this. In addition to the changes in the number of e-folds, 
the shapes of these potentials are such that, as we discuss 
in the next subsection, the tensor tilt is strongly scale 
dependent. Because of the negative value of the scale 
dependent tensor tilt, the gravitational wave background 
amplitude at direct detection scales becomes lower than 
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FIG. 7: The amplitude of the gravitational wave background O.Gwh 2 at 0.1 Hz plotted against the spectral index of the 
gravitational wave background power spectrum at 0.1 Hz. We determine this spectral index using the inflationary flow equations 
with tit{4>) determined at the scalar field value <j) corresponding to 0.1 Hz. In solid, dotted, dot-dashed and dashed curves, we 
show regimes in the J7gw/i 2 — Ut plane at 0.1 Hz populated by power-law, chaotic, symmetry breaking, and hybrid models 
of inflation. In contrast to Figure 1, where there is a large overlap between Monte-Carlo models and analytical models at 
k — 0.002 Mpc -1 , we find that the predictions based on inflationary flow equations depart significantly from predictions of the 
analytical models at k ~ 6.5 x 10 13 . 



the case where tensor tilt is scale independent. 

B. Predictions for Gravitational Wave 
Observations 

With observables at CMB scale defined and models se- 
lected based on the WMAP+SDSS constraints, we now 
extend the discussion to consider the predictions for grav- 
itational wave detection at direct detection frequencies. 
Figure 5 shows the prediction for f^cw^ 2 as a function 
of frequency, where in addition to the power spectrum of 
tensor perturbations generated during inflation, we also 



take into account the evolution through the transfer func- 
tion (See, Section IIC). As can be seen in Figure 5, the 
potentials that have large variations lead to significant 
spectral indices for the tensor power spectrum, while for 
''CMB < 10 -4 the tensor tilt is insignificant and the ap- 
proximation that tit ~ is reasonably accurate. The 
correspondence between large tcmb and large tit impacts 
planning of direct detection experiments. For reference, 
in Figure 5, we also plot the sensitivity curves of three 
of the experimental options that are routinely discussed 
in the literature 0, 0, 0, . The sensitivity curves 
are the same ones shown in Ref. [2(| which were used 
to estimate how well these interferometers can study the 



12 




FIG. 8: The amplitude of the gravitational wave background at 0.1 Hz plotted in terms of the spectral index of the gravitational 
wave background power spectrum at 0.1 Hz. The data points are same as those in the previous figure. Here, we put the 
measurements in the context of experimental studies. The horizontal lines show the minimum amplitude of the gravitational 
wave background detectable with a signal-to-noise ratio of one. The shaded areas represent the regions in the Qcwh 2 — nr 
plane such that these detectors can detect both the amplitude and the tensor spectral index tit at la confidence level. 



gravitational wave background amplitude and the spec- 
tral index. 

Figure 6 shows the relation between the amplitude of 
the gravitational wave background at a frequency of 0.1 
Hz and the tensor-to-scalar ratio at CMB (k = 0.002 
Mpc^ 1 ). Such a relation has been previously discussed 
in Refs. 0, . The data points are distributed such 
that one clearly notices an upper limit on JIqw^ 2 at a 
fixed r, which we determine to be 

S] GW ft 2 < 6.72xl0- 15 r ( P ^mb) \ ( A ™ ) . 
GW V 2 -45 x 10"V V2-74 x 10- 6 J 

This limit agrees with the result of Ref. 17] when con- 
verting our scalar power spectrum normalization and the 



transfer function to their values. This upper limit comes 
from the fact that under slow-roll inflation, hit has a 
negative value if it is not zero. In Fig. 6, for comparison, 
we also produce the prediction under the power-law po- 
tential for inflation. At low values of the tensor-to-scalar 
ratio at CMB, there is a tight correlation between f^cw^ 2 
and f cmb j which falls along the power-law expectation. 

With rcMB < 0.38 at the 95% confidence level from 
current observations, the derived upper limit on the grav- 
itational wave background at 0.1 Hz, under standard 
slow-roll inflation, is such that ficw^ 2 < 2.55 x 10 -15 . 
Experimental options, such as the BBO-lite [2(j, with a 
limiting sensitivity of 10 -15 in flcy/h 2 , are clearly unfa- 
vored as they will only explore an extremely small param- 
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eter space. Moreover, CMB experiments such as Planck 
will soon explore the tensor-to-scalar ratio down to a 
limit of 0.1, while the next generation CMB polariza- 
tion mission targets to improve this down to a limit of 
0.01. In terms of the amplitude alone, a direct detec- 
tion experiment must have sensitivity below 6 x 10~ 17 
in JIqw^ 2 to be both useful in terms of model selec- 
tion and studies related to inflation. Furthermore, even 
if CMB experiments were to detect a primordial tensor 
component with a tensor-to-scalar ratio above 0.01, an 
experimental option such as the standard-BBO [2|| is 
not preferred as many of our inflationary models are such 
that with increasing tcmb , nr (which is negative) also in- 
creases in magnitude. This results in a large suppression 
of power at direct detection frequencies as shown in Fig. 5 
for UgwA 2 versus frequency for curves with r > 0.01. 

The preference to select a direct detection experiment 
with adequate sensitivity down to a CMB tensor-to- 
scalar ratio of 10 -4 is simply a criteria based on slow- 
roll inflationary models. It could be that, under non- 
standard hypotheses, one can generate blue-tilted mod- 
els for tensor fluctuations such that one generates a large 
gravitational wave background at frequencies around 0.1 
Hz. We refer the reader to Ref. 0| for a discussion of 
such possibilities. We, however, find less motivation to 
design a target for a direct detection experiment based 
on non-standard criteria alone since existing cosmological 
data rules out most of these cases. 



C. Experimental Detection 

Figures 7 and 8 highlight the gravitational wave am- 
plitude, in Dgw^ 2 as a function of the tensor tilt at a 
frequency of 0.1 Hz. In addition to the Monte-Carlo 
models generated from the inflationary flow equations, 
we also show the model predictions obtained from ana- 
lytical models. These calculations reproduce results in 
Ref. Figure 8 plots the same information as in Fig- 
ure 7, except that in this plot we have put the model dis- 
tribution in the context of observational measurements. 
The horizontal lines in Figure 8 show the amplitude of the 
gravitational wave background that lead to a signal-to- 
noise ratio (SNR) of one with each of these experiments. 
In terms of the amplitude, BBO, BBO-grand and Ulti- 
mate DECIGO are expected to detect gravitational wave 
backgrounds with amplitudes HqwIi 2 = X x 10 -18 where 
X = 10, 1,0.01 respectively, consistent with previous dis- 
cussions regarding these experiments |l6l l26j . 

In addition to the amplitude of the gravitational wave 
background, it is also useful to consider the possibility 
that these detectors can be used for a reliable measure- 
ment of the tensor spectral index. Following the Fisher 
matrix calculations in Ref. |26j . we determine that for 
a one-sigma detection of the tensor spectral index in 
addition to the gravitational wave background ampli- 
tude, it is required that one roughly satisfies the relation 
\n t \(flGwh 2 /6 x 10~ 18 ) = X. As shown by the shaded 



regions in Figure 8, the requirement that laser interferom- 
eters detect both the amplitude of the gravitational wave 
background and the tensor spectral index leads to a sig- 
nificant reduction in the model parameter space probed 
by these detectors. Models that have both large f^ow^ 2 
and tit may be studied reliably. But even experiments 
such as DECIGO will not be able to detect values of 
\tit\ less than 10 -3 despite their ability to detect the 
gravitational wave amplitudes for nearly all Monte-Carlo 
models. 



D. Testing Inflation: CMB vs. Direct Detection 

While the above discussion considered the extent to 
which the gravitational wave background can be mea- 
sured with planned interferometers such as BBO and 
DECIGO, we now address how CMB information can 
be put to use in determining the prospects for direct de- 
tection measurements. This is motivated by the fact that 
ongoing and upcoming CMB polarization measurements 
will improve the existing constraints from WMAP on the 
CMB tensor-to-scalar ratio and the scalar tilt and run- 
ning of the scalar tilt. With predictions from CMB data, 
one can also optimize the sensitivity requirements for a 
detection of the gravitational wave background. Here, 
we address the extent to which CMB observations can 
be extended to study the gravitational wave background 
at a scale 16 orders of magnitude larger than CMB scale. 

As a first approximation, one can make use of a simple 
power-law extrapolation to establish the expected grav- 
itational wave background at direct detection frequen- 
cies [l4|. To test the accuracy of such an extrapolation, 
we compare the predictions using exact numerical flow 
equations with the predictions calculated using a simple 
power-law extension between CMB and direct detection 
scales. For this comparison, we first write approximate 
forms for the scalar and tensor power spectra as 

Ps{k) » P s (k CMB ) 

( k \ ns (' CCMB ) _1+ 5 Q s(fccMB)ln( X ji s -) 

V &CMB / 

(, \ ii t ():cmb)-1+50t(J;cmb) 
1 " ( 28 ) 

KCMB / 

We can write the gravitational wave background at direct 
detection frequencies as 

fl GW h 2 = A GW P T {k = 0.1Hz) 
P T (k = 0.1Hz) = r . 1Hz Ps(k-0.1Hz). (29) 

We assume that the tensor background satisfies the slow- 
roll consistency relation of ni = — r/8 and that the ten- 
sor power spectrum is not running with aj- = 0. This 
allows us to write an approximate equation for the tensor- 
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FIG. 9: flGwh 2 at 0.1 Hz plotted against the estimated Hgw'' 2 assuming a power-law expansion. The solid line at an angle 

The left panel shows the relation for all data points, while the right 



J power — law • 



represents the case where J1gw/i 2 = [f^Gw/i 2 
panel selects points with 0.99 < n s < 1.01 at CMB scales. The power-law approximation is based on an extrapolation of 
CMB observations using the values rcMB, ns(kcMB), as(kcMB), and the assumption that m = — »*cmb/8 and qt = at all 
scales. This fails to capture the amplitude of the background predicted in our exact numerical calculation. The power-law 
approximation leads to a general overestimate of the gravitational-wave background at direct detection frequencies. While 
some of the scatter in the left panel is due to the large range of n s , by selecting inflationary models with a small range in n s , 
we show that the overestimate is not restricted to either large or small n 3 models at CMB scales. 



to-scalar ratio at direct detection scales as 

nxmz ~ fCMB (30) 

/ l. \ ns(fccMB)-l+|as(fecMB) ln( j +r(fe CMB )/8 



Using this relation and the approximate scalar power- 
spectrum with ns (fccMB ) andas(fccMB) in equation (28), 
one can approximate the gravitational wave background 
amplitude, [^GW^ 2 ]powor-iaw at 0.1 Hz. In Fig. 9, we 
show the comparison between [^cw^ 2 ] power-law and the 
exact amplitude of the gravitational wave background, 
Hqw'' 2 , as calculated through numerical solutions to the 
flow equations. Relative to the exact f^cw^ 2 values that 
are generally below 10~ 15 , [^GW^ 2 ]powcr-iaw ranges as 
high as 10 -10 and over-predicts the gravitational wave 
background at direct detection frequencies. The differ- 
ences are mostly significant for models that involve large 
CMB tensor-to-scalar ratios. These models have a small, 
but non-negligible, scale dependence in the tensor tilt. 

Note that tit = — r/8 and olt — are simple assump- 
tions that do not hold for general inflationary models, 
except in the case of power-law inflation These two 
assumptions also hold well for hybrid models. They do 
not hold, however, for chaotic and symmetry-breaking 
potentials as they often produce significant tensor run- 



ning (see Figure 12). The predictions under the Monte- 
Carlo slow-roll inflation models depart significantly from 
analytical models where olt is mostly zero at direct detec- 
tion scales. This departure is also related to the reason 
why the gravitational wave background amplitude can- 
not be predicted by extrapolating CMB information to 
a scale roughly 16 orders of magnitude higher as shown 
in Figure 9. As highlighted in Figure 7, with non-zero 
ax, chaotic and symmetry-breaking models come closest 
in approximating the slow-roll numerical predictions at 
0.1 Hz, though the agreement is still small. Compared 
to all these analytical models, numerical models lead to 
a substantial scale dependence in the tensor tilt. 

This disagreement is a concern for a simple comparison 
between CMB observations and direct detection experi- 
ments. At CMB scales, the observations are limited to 
TcMB,ns(kcMB) and as(fccMB)- A direct detection ex- 
periment will measure f^cw^ 2 and n/r(0.1Hz). This lim- 
ited information cannot establish a unique inflationary 
potential. Suggestions have been made that the ampli- 
tude of tensor components can be used to establish the 
slope of the inflaton potential V{4>) between CMB and 
direct detection scales. While this is, in principle, possi- 
ble, the constraints will not be unique unless a particular 
class of analytical inflationary models is advocated. Un- 
der the general scenario, where V(4>) can vary such that 
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FIG. 10: Left: The consistency relation between the tensor amplitude r and the tensor tilt tit- Under first order slow roll 
expansion, r = — 8nr- Here, we plot R = r/(— Snr) as a function of tit for measurements at CMB scales. This consistency 
relation is challenging to measure with CMB data alone due to lack of information in the CMB power spectrum for a reliable 
measurement of tit- The spread in R and the tendency for R to dip slightly below one for large nt is a result of the second 
order terms used to calculate r and tit- For —nr > 0.01, one sees a departure from R = 1 due to the fact that for these models, 
tlx > — 2e and the e 2 term is more significant for nr than for r. Right: In a recent study, Smith et al., it was suggested that 
the inflationary consistency relation can be tested by combining CMB observations with direct detection measurements under 
the assumption that nr(fccMB) ~ wt(O.IHz). Here, we plot the tensor spectral index at 0.1 Hz versus the same spectral index 
at CMB showing that nr at 0.1 Hz is larger than the value at CMB scales and the assumption that nr(fccMB) ~ «t (0.1Hz) is 
not consistent with the models predicted under slow-roll inflationary flow equations. In general, using direct detection slope to 
approximate tensor tilt at CMB scales results in a bias of an order of magnitude or more. We explain this difference in terms 
of the scale dependence of the tensor tilt and higher order corrections such as qt. 



the spectral indices are scale dependent, one cannot sim- 
ply connect CMB scales with direct detection frequencies 
through a power-law extrapolation. 

This complication arises due to the large difference in 
the two scales. This difference also complicates simple 
tests on slow-roll inflation. As an example, we revisit 
the proposed test in Smith et al. on the consistency 
relation of r — —8ht, where they assumed nx (^cmb) ~ 
nx(k = 0.1Hz). They made this assumption because 
CMB observations are unlikely to measure the tensor tilt 
with adequate accuracy due to confusion with lensing 
B-modes, while direct detection experiments, especially 
with adequate sensitivity, can make a measurement of 
hit- Using the assumption that fiT(fccMB) ~ nx{k = 
0.1Hz) combined with the tensor-to-scalar ratio at CMB 
scales, Smith et al. advocated a way to test the consis- 
tency relation of the form tcmb = — 8nx(fc = 0.1Hz), 
though resulting constraints for experiments under con- 
sideration were not that impressive. For example, the 
combination of Ultimate-DECIGO and Inflation Probe, 
at best, would constrain this relation to a level of 20% 
if the tensor-to-scalar ratio is between 0.01 and 0.1. The 



test of the consistency relation is affected by the fact that, 
relative to the fractional error on r at any wavenumber, 
one must measure ut with a fractional error that is 8 
times better. 

Beyond these measurement issues, there are also fun- 
damental problems from the theory side. To study 
the extent to which the assumption that titCccmb) ~ 
nrik = 0.1Hz) is accurate, in Fig. 10 (right panel), 
we plot nx(k = 0.1Hz) versus ^(fccMB)- As shown, 
there is a large difference between the two tilts such that 
\nr{k = 0.1Hz)| > |nr(fccMB)|- The consistency relation 
will be violated if one blindly substitutes rir(fccMB) ~ 
nx(k = 0.1Hz) (see, Figure 11), though as shown in the 
left-panel of Figure 10, our predictions are accurate to 
1% level of the consistency relation at CMB with tqmb 
and n T (fc C MB)- 

To show one aspect from which the complications arise, 
in Fig. 12, we plot the tensor tilt at 0.1 Hz as a function 
of the running of the tensor tilt ax, again at 0.1 Hz. For 
comparison, we also plot model predictions for the run- 
ning of the tilt through analytical models, where some 
models, such as the power-law potential, are such that 
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FIG. 11: The consistency relation between the tensor ampli- 
tude r and the tensor tilt ut, R = — r/8nr as a function of 
the GW background amplitude Q,Gvjh 2 at 0.1 Hz under the 
assumption that riT(kcMB) ~ Ut(O.IHz) and using tcmb to 
estimate R. Since ut(O.IHz) > riT(fccMB) by at least an or- 
der of magnitude, R departs from the expected value of one 
by a factor of 10 or more. The figure demonstrates that one 
cannot naively combine CMB tensor-to-scalar ratio with the 
tensor tilt at 0.1 Hz to test the consistency relation. The 
vertical lines mark the limits related to flGwh 2 for three of 
the experiments while the horizontal line R = 1 marks the 
expectation under first-order slow-roll inflation. 



ut = exactly. This large running is a reflection that the 
tensor tilt is largely scale dependent. Simply using the 
tilt at one wavenumber for the tilt at another wavenum- 
ber will lead to incorrect results. Though the running is 
small, the replacement of tilt at either CMB or direct de- 
tection scales with the value of the other is problematic 
since the two scales are largely different. While previ- 
ous studies hailed the large difference in wavenumber be- 
tween CMB and direct detection experiments as a "bless- 
ing" to further study inflation, for arbitrary shaped V(<fi) 
models that are generated numerically, we find that the 
situation is too complex for the simple comparison. The 
large difference in scale requires a careful understanding 
of how parameters behave in between the scales, since 
the potentials are complicated and do not follow a power- 
law behavior of <f>. And since observations are limited to 
two scales, making a potential reconstruction is difficult. 
This, however, mostly applies to the case where one is 
comparing observations to establish information related 
to an arbitrary V{4>). If one assumes a prior shape for 
V(4>), such as through one of the analytical models, then 
the behavior can be predicted analytically and the two 
scales can be matched to the extent that the data allow. 



Instead of combining CMB and direct detection infor- 
mation to test the slow-roll consistency relation, perhaps, 
it is best to improve the measurement of bt(^cmb) with 
CMB polarization data, by improving the sensitivities so 
that confusing lensing B-modes are separated from pri- 
mordial B-modes j35j, and then using CMB data alone 
to study the consistency relation. On the other hand, 
with the feasibility of a direct detection, it is also use- 
ful to investigate if there is way to use information at 
direct detection frequencies for a test of inflation. Un- 
fortunately, due to lack of information on the scalar or 
density perturbations at k ~ 10 13 Mpc -1 it is unlikely 
that information on the tensor power spectrum at di- 
rect detection frequencies alone will be useful to establish 
slow-roll inflation. On the other hand, the complications 
in connecting observations to an underlying model we 
have outlined above may start to improve if information 
on scalar perturbations becomes available at a wavenum- 
ber different from that of CMB but between CMB and 
direct detection scales. Then there would be additional 
data points to interpolate between the two largely sep- 
arated scales. In this respect, in an upcoming paper, 
we will return to such a study by concentrating on the 
power-spectrum of 21-cm background at z ~ 50 to 100 
that provides information on the primordial power spec- 
trum out to k ~ a few tens Mpc -1 . 



IV. SUMMARY 

In addition to density perturbations, inflationary mod- 
els of the early universe generally predict a stochastic 
background of gravitational waves or tensor fluctuations. 
By making use of the inflationary flow approach for single 
field models and with predictions normalized to results 
from cosmic microwave background and large scale struc- 
ture data, we discuss the expected properties of the grav- 
itational wave background from inflation at scales corre- 
sponding to direct detection experiments with laser in- 
terferometers in space. We complement the Monte-Carlo 
calculations based on the inflationary flow equations by 
including predictions expected under several classes of 
analytical inflationary models, including models involv- 
ing hybrid inflation. We find that an improved version of 
Big Bang Observer (BBO-grand) can be used to detect a 
gravitational wave background corresponding to a tensor- 
to-scalar ratio above 10~ 4 at CMB scales, which is two 
orders of magnitude below the published goal of a tensor- 
to-scalar ratio of 10~ 2 by the planned next generation 
CMB polarization space mission by NASA (CMBpol). 

A less sensitive version of BBO (BBO-lite) with a sen- 
sitivity to gravitational waves ficw^ 2 above 10~ 15 is 
unlikely to be useful given that the predictions suggest 
backgrounds well below this level when the tensor-to- 
scalar ratio is below 0.3. Even if the tensor-to-scalar 
ratio were to be above 10 -2 , we suggest that BBO-grand 
will be useful to study inflationary models as standard 
BBO will only allow a marginal detection of the ampli- 
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FIG. 12: The tensor spectral index tit versus the running of the tensor tilt, qt at 0.1 Hz. The data points show the inflationary 
flow predictions, while the shaded regions show the expectations under analytical model with solid, dotted, dot-dashed, and 
dashed lines showing power-law, chaotic, symmetry breaking, and hybrid inflation models, respectively. In certain analytical 
models, such as power-law, olt = over all scales. While the left panel shows olt in a linear scale, the right panel concentrates 
on the low olt region with values around 10 -5 . 



tude while leaving the spectral index unconstrained. The 
polarized foregrounds, such as dust and synchrotron, are 
expected to limit CMB studies to a tensor-to-scalar a 
few times 10 -4 , which is well above the ultimate limit 
of cosmic shear. Unless direct detection experiments are 
also affected by foreground gravitational wave sources, 
laser interferometers may allow us to expand the search 
for gravitational waves down to a tensor-to-scalar ratio 
of 10" 6 . 

Based on our models, we find that simple power-law 
approximations arc not adequate to extrapolate either 
information at CMB scales to direct detection experi- 
ments or vice-versa. The differences are largely under- 
stood through scale dependence of the tensor tilt pro- 
duced in the slow-roll inflationary models through flow 
equations. We also show that the simple assumption of 
«t(^cmb) ~ nr(O.lHz) cannot be used to establish in- 
flationary consistency relations by combining measure- 
ments at CMB and direct detection scales. We advocate 



additional measurements of the primordial scalar power 
spectrum at large wavenumbers beyond CMB as a way 
to improve model selection and predictions of inflation. 
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